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Abstract 

We explore the use of the method of Maximum Entropy (ME) as a 
technique to generate approximations. In a first use of the ME method 
the "exact" canonical probability distribution of a fluid is approximated 
by that of a fluid of hard spheres; ME is used to select an optimal value of 
the hard-sphere diameter. These results coincide with the results obtained 
using the Bogoliuvob variational method. A second more complete use of 
the ME method leads to a better descritption of the soft-core nature of 
the interatomic potential in terms of a statistical mixture of distributions 
corresponding to hard spheres of different diameters. As an example, the 
radial distribution function for a Lennard- Jones fluid (Argon) is compared 
with results from molecular dynamics simulations. There is a considerable 
improvement over the results obtained from the Bogoliuvob principle. 

1 Introduction 

The method of Maximum Entropy (ME) is designed to solve the general problem 
of updating from a prior probability distribution (which often happens to be a 
uniform distribution) to a posterior distribution when new information in the 
form of constraints becomes available. Of all the distributions satisfying the 
constraints the preferred posterior is that which is closest to the prior in the 
sense that it represents the least change of beliefs PQ. 

This suggests that ME can be used to tackle a different kind of problem. 
Even if a distribution is known and it accurately reflects our beliefs it may 
still turn out to be too complicated to be useful in practice. We may need 
to find a more tractable approximation. The idea is to identify a family of 
tractable trial distributions and then select that member which is closest to the 
"exact" distribution by maximizing the appropriate relative entropy. Notice 
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that this is not a problem that can be tackled using the more restricted version 
of the maximum entropy method usually known as MaxEnt. The reason is that 
MaxEnt only allows updating from an underlying physical measure and not 
from a general prior distribution. 

The purpose of this paper is to develop the ME method as a technique to 
generate approximations. In previous work [2] we have shown that a number 
of well-known variational approximation schemes (the Bogoliubov variational 
method, the Kohn-Sham-Hohcnberg density functional formalism, and even the 
variational method of quantum mechanics of which the Hartree-Fock method is 
an important special case) can be derived as special cases of the ME method. In 
2 the ME method was used to develop a mean-field approximation for classical 
fluids in which the tractable family of distributions is obtained by replacing the 
interactions between molecules by an effective external mean field. The effects 
of the long-range attractions are described well but the short-range repulsions 
are badly misrepresented. 

In section 2 the treatment of short-range repulsions is considerably improved 
by selecting as our tractable model a fluid of hard spheres [3]- [5]. The ME 
method is then used (section 3) to select the optimal value of the hard-sphere 
diameter. This is equivalent to applying the Bogoliuvob variational principle 
and reproduces the results obtained by Mansoori et al. 

An alternative approach to the study of fluids is through thermodynamic per- 
turbation theory which yields results formally similar to the variational method. 
In such schemes the actual interatomic potential u is replaced by uq + 8u, where 
uo represents the strong short-range repulsion and 5u is a long-range attraction 
treated as a perturbation. Remarkably, the structure of the fluid for liquid den- 
sities is dominated by the repulsive interactions. The effects of the attraction 5u 
are averaged over many molecules and do not appreciably affect the correlations 
among molecules. As a result, the first-order perturbation theory is quite ac- 
curate. At lower densities, however, higher-order corrections must be included. 
Since Uq is not itself a tractable potential the usual approach is to replace it 
by a hard-sphere potential u^. Several suggestions of how to separate u into 
uq and Su, of where and how to replace uq by Uh s , and how to choose the best 
hard-sphere diameter have been proposed. The most successful are the theory 
of Barker and Henderson |3] and, particularly, the theory of Weeks, Chandler 
and Andersen (WCA) jjj which succeeds in using the hard-sphere Uhs while 
effectively representing the effects of the soft-core potential uq. For a recent 
discussion of some of the strengths and limitations of the perturbative approach 
see 0. 

A clear advantage of the variational and the ME methods is that the im- 
portant and yet somewhat ad hoc nature of the separation of u into Uq and 5u, 
and of the choice of a hard-sphere diameter are eliminated. On the other hand, 
the variational approach fails to take the softness of the repulsive core into ac- 
count, and this leads, in the end, to results that are inferior to the perturbative 
approaches particularly at high temperatures. 

The traditional variational approach allows one to select a single optimal 
diameter; all non-optimal values are ruled out. But, as discussed in P3E], the 
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ME method allows one to proceed further and quantify the extent to which 
non-optimal values should contribute. In this more complete use of the ME 
method, presented in section 4, the "exact" probability distribution of the fluid 
is approximated not by that of a gas of hard spheres with the optimal diameter 
but by a statistical mixture of distributions corresponding to hard spheres of 
different diameters. This is a rather simple and elegant way to take proper 
account of the fact that the actual atoms are not hard spheres but have a soft 
core. The full ME analysis leads to significant improvements over the variational 
method. 

When faced with the difficulty of dealing a system described by an intractable 
Hamiltonian, the traditional approach has been to consider a similar albeit ide- 
alized system described by a simpler more tractable Hamiltonian. The approach 
we have followed here departs from this tradition: our goal is not to identify an 
approximately similar system but rather to identify an approximately similar 
probability distribution. The end result of the ME approach is a probability 
distribution which is a sum over distributions corresponding to different hard- 
sphere diameters. While each term in the sum is of a form that can be associated 
to a real hard-sphere gas, the sum itself is not of the form exp —/3H, and cannot 
be interpreted as describing any physical system. 

In section 5 we test our method by comparing its predictions for Argon 
gas with the numerical molecular dynamics simulation data obtained by Verlet 
|10| . We find that the ME predictions for thermodynamic variables and for the 
radial distribution function are considerable improvements over the Bogoliuvob 
variational result, and are comparable to the perturbative results. Finally, our 
conclusions and some remarks on further improvements are presented in section 



2 Approximation by hard spheres 

We consider a simple fluid composed of N single atom molecules described by 
the Hamiltonian 



where qn — {Pi,fij i = lj --jN} and the many-body interactions are approxi- 
mated by a pair interaction, u(rij) where = \n — rj\. The probability that 
the positions and momenta of the molecules lie within the phase space volume 
dqN is given by canonical distribution, and 
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P(gjv) dq N 



L e -PH{1N) dq 



(2) 



where 



dq N 



Nlh 3N 
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The difficulty, of course, is that this distribution is easy to write down but 
very difficult to use. We must replace P by a more tractable distribution. To 
proceed one must recognize the two features of the interaction potential u that 
seem to be relevant for explaining a wide variety of fluid properties; they are 
the strong short-distance repulsion and the weaker long-distance repulsion. 

To account for the short-distance repulsion we consider the probability dis- 
tribution corresponding to a gas of N hard spheres of diameter r d - The Hamil- 
tonian is 

N 2 N 

H hs (q N \r d ) = ^ — + Uhs with U hs = S ^u hs (r i:j \r d ) , (4) 

i=l i>j 



where 



/ I \ / for r > r d , . 

I oo tor r < r d 

The corresponding probability distribution is 

P hs {qN\r d ) = ^e-l 3H ^^K (6) 
The partition function and the free energy Fh s (T, V,N\r d ) are 

Z hs = J dq N e H»*».(w) d =l f e-^mv.ivlr,) _ (7) 

Two objections that can be raised to the choice of Phs are, first, that it does 
not take the long-range interactions into account; this is a point to which we will 
return later. Second, and this is a more serious problem, it is not clear that Phs 
is a tractable distribution at all. Indeed, the exact solution to the problem of ./V 
hard spheres is not known. However, there exist analytical approximations that 
are in reasonably good agreement with numerical simulations. We will therefore 
assume that for all practical purposes Ph 8 is a tractable distribution. 

For hard spheres the radial correlation function can be calculated within 
the approximation of Percus and Yevick [111 IT2] or, alternatively, from the 
scaled-particle theory |13|. The equation of state can then be computed in 
two alternative ways, either from the so-called "pressure" equation, or from the 
"compressibility" equation but the two results do not agree. It has been found 
that better agreement with simulations and with virial coefficients is obtained 
taking an average of the two results with weights 1/3 and 2/3 respectively. The 
result is the Carnahan and Starling equation of state, 0-0 



PV 



Nk B T 



= ±+v+iLf£ t (8 ) 

hs (1 - V) 



where 77 = ^npr^ with p = N/V. The free energy, derived by integrating the 
equation of state, is 



F hs (T,V,N\r d )=Nk B T 



, . 477 - 3?7 2 



(9) 
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where At = (2-Kh 2 /raksT) 1 / 2 , and the entropy is 



^•H^L = ^ + ^ « io) 

It must be remembered that these expressions are not exact. They are quite 
reasonable approximations for all densities up to almost crystalline densities 
(about r\ m 0.5). However, they fail to predict the face-centered-cubic phase 
when r\ is in the range from 0.5 up the close-packing value of rj « 0.74. 



3 ME determination of the best hard-sphere di- 
ameter 

Next we address the question of which is the best P; ls . Which is the best 
hard-sphere diameter? According to the ME method Q], the preferred trial 
Phs{<lN\fd) is that which is closest to the "exact" P{q^)- It is found by maxi- 
mizing the relative entropy, 

S[P hs \P] = - J dq N P hs (q N \r d ) log ^p^y < 0. (11) 
Substituting eqs. © and ® we obtain 

S [P hs \P] = p[F - F hs - (U - U hs ) hs ] < , (12) 

where (• • • )h s is the average computed over the hard-sphere distribution Phs{QN\fd)- 
Next rewrite ea. lfl2")l as 

F < Fjj = F hs + (U- U hs ) hs , (13) 

and use 

(17 - U ht ) h . = \j d 3 rd 3 r' [u(r - r') - u hs (r - r'\r d )] n™(r, r') , (14) 

where n^(r, r') = (n^ (r,r'))h s and 

f l ( 2 )(r l r') d = S(r — r,;) 5(r' — Tj) — n(r)h(r') — n(r)5(r — r') , (15) 

is the two-particle density distribution. But Uh s {r — r'\r d ) = for |r — r'\ > r d 

(2) 

while nu!{r,r ) = for |r — r | < r d , therefore 

Fu = F hs + (U) hs with (U} hs = ^NpJd 3 ru(r)g hs (r\r d ), (16) 

where we have assumed that the fluid is homogeneous, rv h ' (r, r ) = n\ s (\r — r |), 
and introduced the radial distribution function gh s defined by 

n%(r) d ^p 2 g ha (r\r d ). (17) 
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Maximizing S [Ph s \P] is equivalent to minimizing Fjj over all diameters r d . 
Thus, the variational approximation to the free energy is 



F(T,V,N)RtFu{T,V,N\r m ) 



del 



mmFu(T,V,N\r d ) 



(18) 



where r m is the optimal diameter. Notice that the approximation does not 
consist of merely replacing the free energy F by a hard-sphere free energy Fh s 
which does not include the effects of long range attraction; F is approximated 
by Fjj(r m ) which includes attraction effects through the second term in ea. H16|) . 
This addresses the first of the two objections mentioned earlier. Indeed, the 
real fluid with interactions given by u is not being replaced by a hard-sphere 
fluid with interactions given by uts] it is just the probability distribution that 
is being replaced in this way. The internal energy is approximated by (H)h s = 
%Nk B T + (U) hs and not by (H hs ) hs = %Nk B T. 

To complete the solution one needs the radial distribution function gh s ( r \ r d)- 
A numerical solution by Throop and Bcarman is tabulated in |14| . Alternatively 
it can be written in terms of the Laplace transform of rg^sfr \r d ) |12| . 



G(s) 



dx xg hs (xr d \r d )e~ 



sL{s) 



,o l2 V [L(s) + S(s)e s 

where x is a dimensionless variable x = r/r d , 



L(s) = \2r] 



s + (1 + 2?/) 



(19) 



(20) 



and 



S(s) = (1 - r)f s 3 + 6»? (1 - v) s 2 + 18?7 2 s - 12r? (1 + 2rj) . (21) 

To make use of this solution it is convenient to write {U)h s in terms of V(s), 
the inverse Laplace transform of rit(r), 



which gives 



xu(xr d ) = ds V(s)e 
Jo 



(U) hs = 12N V / ds V(s)G(s) 



For a Lennard- Jones potential, 



u(r) = 4e 



(22) 



(23) 



(24) 



we have 



V(s) = Ae 



. 12 10 

a \ s u 



Uj 10! \r d J 4! 



> 6 4 

a \ s 



(25) 
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Finally, it remains to minimize Fjj in eci. (|16[) . This is done numerically; an 
explicit example for Argon is calculated in section 5. 

As mentioned in the Introduction a problem with the approach outlined 
above is that it fails to take the softness of the repulsive core into account. This 
flaw is manifested in a less satisfactory prediction of thermodynamic variables 
at high temperatures, and also, as will be shown by the numerical calculations 
in section 4, in a poor prediction of the radial distribution function. However, as 
we shall see in the following sections, the ME method does not stop at the mere 
selection of the best diameter. A more complete ME analysis offers significant 
improvements over the variational method. 



4 A more complete ME analysis 

The ME method as pursued in the last section has led us to determine an optimal 
value of the hard-sphere diameter. Next we ask to what extent do we believe 
that the correct selection should be r d = r m rather than r d = r m + Sr. To what 
extent are values 5r ^ ruled out by the ME method? ^ |U] This question 
is an inquiry about the probability of r d , Pd{i"d)- Thus, we are uncertain not 
just about qN given r d , but also about the right r d and what we actually seek 
is the joint probability of q^ and r d , Pj{qN, r d)- Once this joint distribution is 
obtained our best assessment of the distribution of qjq should be given by the 
marginal over r d , 

Phs(qN) == J dr d Pj(qN,r d ) — J ' dr d Pd{rd)Phs{qN\r d ) (26) 

By recognizing that diameters other than r m are not ruled out and that a 
more honest representation is an average over all hard-sphere diameters we are 
effectively replacing the hard spheres by a soft-core potential. 

Next, we make a second use of the ME method to obtain Pj(qNii" d ). Of 
all distributions within the family Pj{qN ,Td) — P d {r d )Phs(qN\r d ) we seek that 
which is closest to the prior m(qN,r d ) over the space {qN,r d }. But what is 
m(qN,r d )7 We want the marginal distribution Phs(qN), eq.(|2EJ), as close as 
possible to the "exact" distribution P(qN) in eq.©. This is achieved if we 
set m(qN,rd) — P(qN)^(fd)- Furthermore, we can assign [i(rd) noting that 
rd is not a physical quantity, its only meaning derives from being a parameter 
appearing in the distribution Phs{QN\ r d)- Then there is a unique natural distance 
d£ 2 = r y(rd)drj l in the space of r^s which is given by the Fisher-Rao metric, 

( \ vi I , f d\ogP hs {q N \r d )\ 2 , , 

l[ r d) = / dq N P hs (q N \r d ) I ^- I . (27) 

Therefore n{rd) — 7 1 ^ 2 (f<i)- Therefore the joint Pj(qN,Td), or equivalently the 
diameter distribution Pd(rd), is determined by maximizing the entropy 

a [P d ] = - I ' dq N dr d Pd(r d )P hs (q N \r d ) log ?^fhf^M , ( 28 ) 
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For sake of simplicity, one can reorginize this entropy to 

a [Pd] = -J dr d P d (r d ) log + J dr d P d (r d )S [P hs \P] , (29) 

where <S [-P/i S |P] is given eqs.{12| an d (fTHfl . 

S [P hs \P] =P{F-F V ) = P(F- F hs - (U) hs ) . (30) 

Maximizing a [P d ] over variations 5P d subject to a normalization constraint 
gives, 

p S[P hs \P] p-pFv 

Pd(r d )dr d = e — 7 1/2 (U) dr d = — — ^ 2 (r d ) dr d , (31) 

where the partition functions Q and Cu are given by 

C = eP F C v with Ca = / dr d 7 1 / 2 (r d ) e^ F " . (32) 



The remaining problem in the above equations is the calculation of the Fisher- 
Rao measure 7 1 / 2 and this is conveniently done by considering the entropy of 
Phs(qN \r' d ) relative to P hs (qN \r d ), 

S[P hs (q N \r' d ) \P hs (q N \r d )} = - [ dq N P hs (q N \r' d ) log ^ ^ \' d) . (33) 

J Phs{qN\r d ) 

A straightforward differentiation shows that 

d 2 S[P hs {q N \r' d )\P hs {q N \r d )] 



dr'j 



- l{r d ) ■ (34) 

r' d =r d 



Substituting the distributions Ph s {qN \r' d ) and Phs(qN \r d ) into ea. (|33)l gives 



S [Phs(qN\r' d )\P hs (qN M] = (3 [F hs (r d ) - F hs (r' d ) - (U hs (r rf )) r , + (U hs (r' d )) r , d 

(35) 

where (• • -) r / is the average over Ph s {qN\r' d ) ■ As we argued above ea. l(TB|l the 
expectation of the potential energy (Uhs { r ' d ))r' d vanishes because the intermolec- 
ular potential u{r \r' d ) is zero in the region r > r' d where the radial distribution 
function gh s {r \r' d ) is not zero. A similar argument shows that (Uhs (f d )}r d = 
when r' d > r d . However, when r' d < r d the expectation (Uhs (^d))r^ diverges, S 
is not defined and ea. (|34|l is not applicable. We can argue our way out of this 
quandary by pointing out that the divergence is a consequence of the unphys- 
ical idealization involved in taking a hard-sphere potential seriously. For more 
realistic continuous potentials the distance between r' d = r d + dr d and r d is the 
same as the distance between r' d = r d — dr d and r d . We can then always choose 
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r' d > rd and define 7(r<z) in ea. llMI) as the limit r' d = rd + + . Then, using eq.@, 
we have 



l{r d ) = (5 



d 2 F hs (r' d ) 



di 



--r d +0+ (1 - V) 



To summarize, the distribution of diameters Pd{rd) is given by ea. l|31|) with 
Fjj given by eqs. ((TBI ESJ and 7 given by Our best approximation 

to the "exact" P{qn) is the Phs(qN) given in ea. l(^|) . The corresponding best 
approximation to the radial distribution function is 

9hs{r) = J dr d Pd(rd)ghs(r \r d ) ■ (37) 

When intermolecular interactions are given by two-body potentials the equa- 
tion of state, the free energy, the internal energy, and other thermodynamic 
functions can be expressed in terms of the radial distribution function. There- 
fore an improved g(r) leads to improved estimates for all other quantities. 

However, there is a problem. For large N the distribution Pd(rd) ~ exp —(5Fu 
is very sharply peaked about the maximum attained at the optimal diameter 
r m because Fjj is an extensive quantity Fjj cx N. This result must be inter- 
preted with care: when choosing a single optimal diameter for a macroscopic 
fluid sample we find that the ME confers overwhelming probability to the opti- 
mal value. This is not surprising. The same thing happens when we calculate 
the global temperature or density of a macroscopic sample: the ME method 
predicts that fluctuations about the expected value are utterly negligible. And 
yet fluctuations can be important. For example, for smaller fluid samples, or 
when we consider the local behavior of the fluid, fluctuations are not merely 
observable but can be large. Local fluctuations can be appreciable while global 
fluctuations remain negligible. 

The question then, is whether these local fluctuations are relevant to the 
particular quantities we want to calculate. We argue that they are. The radial 
distribution function g(r) is the crucial quantity from which all other thermo- 
dynamic variables are computed. But from its very definition - g(r) is the 
probability that given one atom at a given fixed point, another atom will be 
found at a distance r - it is clear that g(r) refers to purely local behavior and 
should be affected by local fluctuations. To the extent that the optimal r m 
depends on temperature and density we expect that local temperature and/or 
density fluctuations would induce local diameter fluctuations as well. 

It is important to note that unlike density fluctuations, the local diameter 
fluctuations are not real. They cannot be: the hard spheres do not exist except 
in our minds. The diameter fluctuations are merely a representation of our 
uncertainty about which r d to choose. 

The extended analysis in this section does not yet allow us to pursue the 
question of local fluctuations in a satisfactory manner. For the purpose of this 
paper, however, we can quickly estimate the effects of local fluctuations by 
assuming that the effective number of atoms N e g that are locally relevant to 
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the calculation of g(r) is smaller than TV. We deal with an effectively smaller 
fluid sample. The actual calculation of N e g will be pursued elsewhere; it is not 
difficult to see that the ME method itself still applies all that is needed is a 
broader family of trial distributions. 

Next we apply the extended ME method developed above to fluid Argon. 

5 An example: Argon 

One of the difficulties in testing theories about fluids against experimental data 
is that it is not easy to see whether discrepancies are to be blamed to a faulty 
approximation or to a wrong intermolecular potential. This is why it is theories 
are normally tested against molecular dynamics numerical simulations where 
there is control over the intermolecular potential. In this section we compare 
ME results against simulation results 10 for fluid of monatomic molecules in- 
teracting through a Lennard- Jones potential H24|) . The parameters e and a (the 
depth of the well and the radius of the repulsive core, u(a) = 0) are chosen to 
model Argon: e = 10~ 2 eV and a = 3.405 A. 

Figure[I] shows the free energy Fjj /NksT as a function of hard-sphere diam- 
eter Td for Argon at a density of pa 3 — 0.65 for different temperatures, and fig'El 
shows Fu/NksT as a function of r& for several densities at fixed T = 107.2 K. 
Since the critical point for Argon is T c = 150.69 K and density per 3 = 0.33 
all these curves lie well within the liquid phase. The increase of Fu/NksT 
for low values of r d is due to short range repulsion described by Fhs/NksT 
while the increase for large r d is due to the long range attraction described by 
(U) hs /Nk B T. 

The best is that which minimizes Fjj and depends both on temperature 
and density. The best diameter decreases as the temperature increases: atoms 
with higher energy can penetrate deeper into the repulsive core. The dependence 
with density is less pronounced. 

Next we study the distribution of diameters, ea. (f3*TJ) . As discussed earlier 
for large N the distribution Pd{rd) ~ exp — (3Fu is sharply peaked about the 
optimal diameter r m . But we argued that the effective number of particles that 
are relevant to the local behavior is smaller N e g. In Fig|3|we plot the distribution 
Pd{fd) for different temperatures, for a fixed fluid density of per 3 = 0.65, and 
for an arbitrarily chosen fixed N g — 13500. Notice again that the distribution 
shifts to lower diameters as the temperature increases. Notice also that at 
lower temperatures the distribution becomes narrower. This means that the 
effects of the softness of the repulsive core are less important; that a hard-sphere 
approximation is better at low temperatures 3 . 

Finally in fig. (J2J we show a comparison of the radial distribution functions 
computed in three ways. The dotted line is ghs(r\r m ) for the hard-sphere fluid 
with optimal diameter r m . This curve, calculated from the numerical data 
tabulated in ^3] , is also the result of the variational method and coincides with 
the ME result for a macroscopically large N c s — N . The dashed line is Ver let's 
molecular dynamics simulation The continuous line is the averaged g~hs{r) of 
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Figure 1: The free energy Fjj as a function of hard-sphere diameter for 
Argon at a density of per 3 = 0.65 for different temperatures. The best r<j is that 
which minimizes Fjj . 




2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 



r d ( A ) 



Figure 2: The free energy Fjj as a function of hard-sphere diameter r<j for Argon 
at T — 107.82 K for different densities. The best is that which minimizes Fu- 
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Figure 3: The distribution of hard-sphere diameters for Argon for several 
temperatures at density pa 3 — 0.65. The size of fluid was arbitrarily chosen as 
N eS = 13500. 

2.8- Improved ME method with effective N=19171 

2 6 J PY hard sphere solution with r =3.4149A 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 



Figure 4: The radial distribution function for (a) the hard-sphere fluid with 
optimal diameter r m ; (b) Verlet's molecular dynamics simulation; and (c) the 
improved ME analysis, for Argon at density per 3 = 0.65, temperature T = 
107.82 K, and effective particle number N cS = 19171. 
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the more extended ME analysis. The calculation is for a density per 3 = 0.65, 
temperature T — 107.82 K, and the value of effective particle number N e g = 
19171 was adjusted to achieve the best fit. The agreement between the ME 
curve and Ver let's data is remarkably good. The vast improvement over the 
simpler variational method calculation is clear. 

One might be tempted to dismiss this achievement as due to the adjustment 
of the parameter N e g but this is not quite fair: there is a single free parameter 
and the functional form of the whole curve ghs{r) in ea. (|37[) is correctly repro- 
duced. We should point out, for example, that the Fisher-Rao measure 7 1 / 2 (r^) 
plays a crucial role. Omitting the factor 7 1 / 2 (r^) from the distribution Pd{rd) 
would have led to a much less satisfactory fit regardless of the choice of N e g. 

6 Conclusion 

The main goal of this paper has been to show that the ME method can be used 
to generate approximations in a way that generalizes the Bogoliuvob variational 
principle. This addresses a range of applications that lie beyond the scope of 
the traditional MaxEnt. 

We showed that rather than approximating the real fluid by a fictitious one, 
it is better to approximate the "exact" probability distribution by a statistical 
mixture of distributions corresponding to hard spheres of different diameters 
which allows a better description of the soft core of the short-range repulsive 
potential. 

The results achieved in this paper represent a step in the right direction but 
they are by no means final. When we say that the averaged PhsiQN) is the 
"best" approximation to P(qpj) we do not mean that one cannot do better. In 
fact, further improvements are always achievable by choosing a broader family 
of trial distributions. In the study of fluids we saw that these improvements 
are not just possible, they are necessary. We argued that the next important 
improvement of the ME method as a calculational tool for fluids should be in 
the direction of developing a theory of local fluctuations. This would lead to 
a systematic method for the determination of the effective number of particles 
N e g that are locally relevant. 

We conclude, therefore, that the ME approach to the theory of fluids is 
already a significant improvement over the Bogoliuvob variational approach. 
With only moderate further developments we can realistically expect it to sur- 
pass even the best perturbative methods developed to date. 
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